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Abstract 

j^ ■ As a first step to understanding the role of molecular or chemical polydispersity in self-assembly, 

O ■ we put forward a coarse-grained model that describes the spontaneous formation of quasi-linear 

• \ polymers in solutions containing two self-assembling species. Our theoretical framework is based 

Cd I on a two-component self-assembled Ising model in which the bidispersity is parameterized in terms 

Vh ■ of the strengths of the binding free energies that depend on the monomer species involved in the 

n^ • pairing interaction. Depending upon the relative values of the binding free energies involved, differ- 

C . ent morphologies of assemblies that include both components are formed, exhibiting paramagnetic- 

Q I , ferromagnetic- or anti ferromagnetic-like order, i.e., random, blocky or alternating ordering of the 

I two components in the assemblies. Analyzing the model for the case of ferromagnetic ordering, 

CN ■ which is of most practical interest, we find that the transition from conditions of minimal assembly 

JT ■ to those characterized by strong polymerization can be described by a critical concentration that 

l/^ . depends on the concentration ratio of the two species. Interestingly, the distribution of monomers 
in the assemblies is different from that in the original distribution, i.e., the ratio of the concen- 
trations of the two components put into the system. The monomers with a smaller binding free 

f~^ ■ energy are more abundant in short assemblies and monomers with a larger binding affinity are 

^^ ■ more abundant in longer assemblies. Under certain conditions the two components congregate 

. . . into separate supramolecular polymeric species and in that sense phase separate. We find strong 

. ^ I deviations from the expected growth law for supramolecular polymers even for modest amounts of 

^\ \ a second component, provided it is chemically sufficiently distinct from the main one. 
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I. INTRODUCTION 

Self-assembly processes play a key role in the construction of biological structures and ma- 
terials including the cell skeleton, viruses, bone, protein complexes and amyloid fibrils [l|-l6|. 
These natural supramolecular structures have inspired scientists to exploit supramolecular 
self-assembly principles as a tool in order to design and build bio-mimetic molecular struc- 
tures from synthetic molecular compounds for purposes such as drug delivery and biomedical 
diagnostic technologies [7H9|. These novel materials have found applications in nanotechnol- 
ogy, medicine, including dental, cosmetic surgery and orthopedic applications, and personal 
care products [iCHlG 



Application of bio-inspired materials requires mass production of their molecular building 
blocks with efficient methods, which are not necessarily as accurate as the ones employed 
in research laboratories, or, for that matter, biology. Indeed, industrially produced self- 
assembling molecular units tend not to be very monodisperse, i.e., consist only of a single 
compound, but often consist of a large number of similar molecules with varying size, charge, 
chemical composition, and so on. Molecules not chemically identical to the target molecule 
are sometimes called impurities, but the whole collection may also be seen in some general- 



ized sense as a polydisperse one [17|, Il8 



Taking as an example beta-sheet forming peptides, polydispersity of this kind may lead 
to inter-molecular binding affinities that vary as a function of the polydispersity attribute of 
interest. This may in this case include amino acid sequence, the number of amino acids that 
make up a peptide, small molecule reaction products that are able to hydrogen bond to the 
peptides and so on [12]. All of this may have a large impact on the solution properties of 
the assemblies and on the structure of the assemblies themselves. Therefore, understanding 
the role of polydispersity on the nature of order- disorder transitions and morphologies of 
spontaneously formed supramolecular structures is very important from both a fundamental 
and an applied scientific point of view. 

Despite the fact that self-assembling molecular blocks are always to some degree chem- 
ically polydisperse, this issue has, as far as we are aware, received little attention in the 
literature. The influence of polydispersity has been studied theoretically in the context of 
self-assembling block copolymers |19|-|22|. and theoretically and experimentally in that of 



chiral ampliflcation in supramolecular polymers |23H30| . In the latter, the net observed 
chirality of a solution is studied by varying the composition in mixtures of enantiomers of 
self-assembling compounds and in mixtures of achiral and chiral species that co-assemble 
J24J . l3l| . Because studies like these focus entirely on the net macroscopic helicity, very little 
information is available on the structure and composition of the assemblies. 

The aim of the present work is to get insight into the role of polydispersity in the simplest 
case of self-assembly, i.e., that of quasi-linear self-assemblies also called equilibrium polymers 
or EPs. In EPs all the monomers are able to bond reversibly with each other with a binding 
free energy and form quasi-linear self-assemblies of varying length. The growth of these 
assemblies usually takes place in the form of nucleated assembly in which, due to a required 
change of conformation of monomers in the bound state, an activation free energy barrier 
must be overcome in order to form assemblies. As one increases the density of free monomers 
in the system, one observes a transition from a regime of mainly monomeric state to a self- 
assembly dominated regime. The transition from one regime to another can be characterized 



by a critical concentration that depends on both nucleation and binding free energies [32M34 



An example of such systems is given by beta-sheet-forming peptides in which the binding 



between the peptides takes place through hydrogen bonding between the oxygen in the 
backbone of the amino acid residues in one peptide and nitrogen in the backbone of amino 
acids in the other peptide (g, [l2|. The activation free energy in this case originates from 
the transformation of conformation of peptides from a random coil or alpha-helical state 
with a lower free energy to a an extended, rod-like conformation with a higher free energy. 
Nevertheless, peptides joined to a beta-sheet gain a binding free energy that is large enough 
to counterbalance the free energy of the rod-like bound state. 

As a first step towards understanding the influence of polydispersity in linear self- 
assemblies, we focus our theoretical study on a bidisperse self-assembling system, schemat- 
ically illustrated in Fig. 1. To model linear self-assembly in such systems, we propose a 
two-component lattice-gas model in which bidispersity is incorporated through the species- 
dependent free energy parameters describing the binding and nucleation processes ^] . Map- 
ping the problem onto a 1-D Ising model and invoking the standard transfer matrix method, 
we calculate the partition function and explore the ordering of the two components in as- 
semblies of arbitrary length. This allows us to probe the composition assemblies of all sizes 
as function of concentration and interaction strengths between bound monomers. 

We note that, in spirit at least, our approach is similar to that found in other works 
in the wider context of polymer physics, including the helix-coil transition of polypeptides 
and the melting transition of DNA |26l . 1271 . |35|, as well as structural transformations in 



two-component copolymers |28|, where transfer matrix methods have been employed. Fur- 
thermore, lateral ordering in tape-like structures of annealed blocky copolymers has been 



studied by simulations [29|, |30[. The difference with the earlier work lies in the coupling 
between self-assembly, composition and length distribution of the assemblies. 

The binding free energies describing our equilibrium polymerization model are denoted 
bij > 0, measuring the free energy gain (in units of thermal energy) of the bonded interaction 
between two monomers of type i,j = 1,2, hence the positive sign, where 612 = &21 by 
symmetry. For definiteness, we presume that bu > 622, so molecules of species 1 forge 
stronger bonds to each other than species 2 do. There are also two activation free energies 
ai associated with the two species z = 1, 2 describing the (dimensionless) free energy penalty 
associated with conformational changes when monomers are absorbed in assemblies. See 
Fig. 1. 

Our aim is to answer the following two three central questions: 

i) What happens to the mean size of the assemblies if we mix two distinct self-assembling 
species that are able to co-assemble? 

ii) Do the two components actually co-assemble into linear aggregates, or do they form 
chemically pure linear assemblies consisting of one species only? 

iii) What physical principles regulate the composition of the assemblies? 

We find that the relevant quantity determining the arrangement of monomer species in the 
assemblies (the morphology) is an effective "coupling constant" J = |(6ii + 622 — 2612) in 
the parlance of the Ising model, which depends on a linear combination of the binding free 
energies bij. This coupling constant is, of course, not unlike the Flory-Huggins parameter in 
binary lattice fluids |36| Depending on the value of this effective coupling constant J, three 
different morphologies can be envisaged. For J > 0, the two types of monomer are organized 
in linear self-assemblies in a blocky, "ferromagnetic-like" order. The "paramagnetic" J = 
case corresponds to a random distribution of monomers in the assemblies, whereas J < 



leads to "anti- ferromagnetic-like" , alternating ordering of the monomers along the assem- 
blies. In the limit J ^ — 1, the latter would represent co-ordination polymers |37|. The 
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FIG. 1: A schematic of the linearly self-assembling system under study and the morphology of 
the structures formed. The two species of self assembler give rise to two activated states with 
associated free energies oi and a2, accounting for conformational changes necessary for binding, 
and three binding free energy gains bu, 612 = ^21 and 622 



"paramagnetic" (random) case corresponds to a self-averaging system that exactly behaves 
like a monodisperse system and therefore is not of interest in our work. The "ferromagnetic" 
(blocky) case occurs when the binding-free energy between the two distinct species is less 
than the average binding free energy of the two species. This is the case that is relevant to 
those situations we are interested in, i.e., those cases where polydispersity as defined earlier 
plays a role. For instance, for the beta-sheet forming oligopeptides discussed earlier, and 
reported on by Aggeli and collaborators, binding takes place through hydrogen bonding 12 . 
When a short peptide binds to a longer one, the number of hydrogen bonds is equal to that 
formed between two short ones, suggesting that 612 = &22 and J = {bu — 622 )/4 > 0, so this 
would indeed correspond to the "ferromagnetic" (blocky copolymeric) case. 

Although the "ferromagnetic" case J > does obviously not cover all possible circum- 
stances, one can envisage many experimental situations where this condition is met. For 
completeness we present our results for arbitrary values of J throughout the paper although 
we do analyze our model and discuss the consequences of bidispersity only for the ferromag- 
netic case. According to our findings, the bidisperse system behaves in many ways similar 
to the monodisperse one, and exhibits a transition from a regime with minimal assembly to 
one where self-assembly predominates. If sufficiently co-operative, the polymerization tran- 
sition is sharp and is demarcated by a critical concentration or temperature. The precise 
value of the critical concentration or temperature depends on the stoichiometric ratio of the 
concentrations of the two components involved. 

Furthermore, for large asymmetries between the two species, i.e., large J values, and 
for sufficiently low concentrations, a "demixing" region appears where pure self-assemblies. 



composed mainly of one species, coexist. Note that there is no demixing here on macroscopic 
scale, only on the scale of the individual assemblies. Considering the distribution of the 
two species along the assemblies, we find that the distribution of the two species in the 
assemblies differs from that of the parent distribution as the size of the assemblies grows. 
In short assemblies the population of species with less binding affinity is dominant, while in 
longer assemblies monomers with a bigger binding affinity are more abundant. 

The remainder of this paper is organized as follows. In Section 2, we first outline a 
theoretical framework that is based on a two-component self-assembled lattice-gas model. 
In section 3, we show that mapping our model onto an Ising model allows us to explore and 
predict the "phase" behavior of bidisperse monomers forming quasi-linear self-assemblies. In 
Section 4, we explore the self-assembly behavior in the low- and high-concentration limits. 
Sections 5 and 6 are devoted to a comprehensive analysis of our model, focusing in particular 
on a bidisperse system with ferromagnetic-like (blocky copolymeric) ordering. We discuss 
and compare the results for the cases of small and large J. Of particular interest is the 
relative fraction of the two species in the assemblies as their size varies. Finally, we conclude 
our work in section 7, where we summarize our main findings and discuss and compare the 
infiuence of bidispersity in self-assembling systems with other thermodynamical systems. 



II. EQUILIBRIUM STATISTICS OF BIDISPERSE SELF-ASSEMBLING 
MONOMERS 

We consider a model system consisting of two self-assembling species that can form lin- 
ear self-assemblies of any length A^ = l,2,...,oo . Each monomer in a typical linear self- 
assembly can be either of species 1 or 2. The difference between the two species arise either 
from their size, chemical structure and so on. We parameterize this difference between the 
species by invoking effective activation and binding free energies that depend on the type 
of species. Again, let hij > denote the free energy gain of the bonded interaction between 
two monomers of type z, j = 1, 2 in a self-assembly, and Oj > the activation free energy of 
a monomer of type i. In principle, the value of the former depends on the next neighbors 
of the two binding molecules but we ignore this complication here. In our model, bij = bji 
and tti > 0, where for definiteness we suppose bu > ^22- All free energies are scaled to the 
thermal energy kBT = 1. Here, fc^ denotes the Boltzmann constant and T the absolute 
temperature. 

More generally, these effective free energies can originate from different types of interac- 
tions such as electrostatic, hydrogen bonding, hydrophobic and van der Waals interactions. 



For instance, for the pertinent system of Aggeli et al. mentioned earlier [12| , involving pep- 
tides of L amino-acids as monomers, the binding free energy gains bij result from hydrogen 
bonds between the oxygen atoms in the backbone of one peptide and the nitrogen atoms 
in the backbone of the other peptide. This gives rise to 2 hydrogen bonds per residue and 
2L hydrogen bonds in total, suggesting that ba = boL, with bo a proportionality constant. 
Similarly, 612 results from the hydrogen bonding between the backbones of two peptides of 
different length, and the number of hydrogen bonds in this case is equal to twice the number 
of residues of the shorter peptide, so 612 = &22- 

The transformation free energies Cj result from loss of conformational entropy due to a 
change of configuration of peptides from a single, coil-like monomer to a rod-like extended 
and a bound state in the assembly. As a first order approximation, one would expect them 
to depend linearly on the size of oligopeptides, i.e., a^ = aoLi. 



The (diinensionless) grand potential energy of a system consisting of free monomers and 
self-assembled polymers can be written as the sum of an ideal entropy of mixing and the 
contribution of the internal partition functions of all the assemblies of varying size A^, giving 

- = ^p(iV)[ln(p(iV)z/) - 1 - lnZjv(p„6.„a,)], (1) 

N=l 

where V is the volume of the system, Z^ the (semi-grand partition) function of an assembly 
consisting of N monomers. Here, z/ denotes the interaction volume equal roughly to an 
effective volume of a solvent molecule [3^, which we assume to be independent of the 
type of species, and p{N) the number density of self-assemblies of size N. The semi-grand 
partition function Zi<^ counts the number of configurational states of linear assemblies of size 
A^, composed of the two species of monomer and described by their relevant binding and 
transformation free energies, and dimensionless chemical potentials fii {i = 1,2). The latter 
are fixed by the total concentration of the monomers in the solution. Our grand potential 
tacidly assumes the assemblies to be in the dilute limit, i.e., interactions between assemblies 
are presumed to be negligible. 

The equilibrium size distribution minimizes the grand potential Eq. [1] Setting 
6n/6p{N) = yields 

p(N) = v'^Zn{pu bij, ai). (2) 

Therefore, our task of finding the equilibrium size distribution reduces to the calculation 
of the semi-grand partition function of a polymer with A^ degrees of polymerization, whose 
monomers can be either of the two species. To do this, we model a linear self-assembly of 
length A^ as an one-dimensional interacting two-component lattice gas of size A^ each site 
of which is occupied by either of species 1 or 2. Each site numbered / = 1, ..., A^ can be 
identified by the occupation numbers n| = 0, 1 and nl = 1 — nj describing the number of 
monomers of species 1 and 2 at position / along the lattice, respectively. These occupation 
numbers automatically obey the constraint E^;^(?7,J + nl ) = N. 

The mutual exclusivity of occupation of each site by either of the species gives us the 
possibility of mapping the two occupation numbers to a single "spin" variable of the Ising 
model of ferromangetism by a simple transformation 
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(1 + 50, (3) 



nP = lil-S,), (4) 

where Si are the "spin" variables that take ±1 values. As a result, in this representation 
the species 1 can be though of as an "up" spin and the species 2 as a "down" spin. In order 
to calculate the semi-grand partition function Z^, we first express the Gibbs free energy of 
a self-assembly of length A^ in terms of the spin variables. The (dimensionless) internal free 
energy of monomers of type i in free solution is lower by an amount equal to —a^ compared 
to those of bonded ones due to the larger number of available degrees of freedom in the 



unbound state. 




G(l) = -|(l + 50-|(l-50 


(5) 


G{N > 1) = )J[-^(1 + ^.)(1 + ^m) - ^(1 - S,){1 - 5.+1) 


(6) 



i=l 



-^{(1 + ^.)(1 - ^.+1) + (1 - ^.)(1 + ^.+1)}] 



The semi-grand partition function can then be written as 



N 



{Si} 1=1 

The partition function for monomers in the unbound, solution state has the simple form of 
Zi = exp(/ii + fli) + exp(/i2 + 02). 

For all cases A^ > 1, we find, by carrying out the sum over the spin values and simplifying 
the semi-grand partition function by collecting terms of equal order, that it is equivalent to 
the partition function of an Ising chain of size N. Its Hamiltonian reads 



N-l N 

nN>l{{Sl}) = -JY1 SlSl+^ -HJ2^1- ^o(^) + i(^ll - ^22)('5l + Sn), (8) 

;=i 1=1 

with the effective coupling constant J, magnetic field strength H and an energy term eo{N) 
that is an invariant of the spin states, defined as 

J = -{bu+b22-2bu), (9) 

i^ = ^[(611 -M + (-"1-/^2)], (10) 

SoiN) ^ ^^^^(&ii + &22 + 2612) + y (/^l + /^2), (11) 



In the special case of 622 = &12, relevant for the oligopeptides of Aggeli et al. [12[ discussed 
earlier, we obtain J = l{bu - ^22) and eo{N) = {N - l){bu + 3622)/4 + iV(/ii + /i2)/2. 

If we introduce the renormalized chemical potentials /i^ = fii + bu, this gives H = (/xi— /i2) 
and eo(iV) = N{—J + {fi[ + /i2)/2) — b, where b = (611 + 622 + 2&i2)/4 is as before the mean 
binding free energy. We conclude that there are only two relevant energetic parameters, J 
and b. The average binding free energy does not couple to the spin states, is independent 
of A^, and plays a role similar to that of the binding free energy in the monodisperse case. 
The parameter J does couple to the spin states of the lattice and determines in the end 
the composition of an arbitrary assembly. In the next section, we describe how to calculate 
the partition function of our model and other relevant quantities such as fraction of self- 
assemblies by exploiting the mapping onto the Ising model. 



III. MAPPING ONTO THE ISING MODEL 

Mapping our problem onto the one-dimensional Ising model already provides us with 
a direct insight into the general "phase behavior" of the system at hand - here "phase 



behavior" does not refer to phase separation on the macroscopic scale but rather on the 
microscopic scale, that is, between assemblies. As already advertised in section I, depending 
on the value of the coupling constant J, associated with the spin (or occupation) states of 
two neighboring sites, different types of ordering may appear in the assemblies. For J > 
ferromagnetic ordering is favored, implying blocky copolymers, for J = the paramagnetic 
case is favored, meaning random copolymers, and for J < anti-ferromagnetic ordering, 
associated with alternating, copolymeric ordering. See Fig. [T] 

Exploiting the standard method of the transfer matrix |38| , the resulting partition func- 
tion for free boundary conditions (implying no preference for any monomer to sit at the ends 
of the assemblies) takes the form, 

Zn>i = [a;+A^-^ + x_A^-i]exp(£o(iV)), (12) 

(1 + e2^)e2^ ± v/4e2^ + [e^" - lye^J 
^± = ^^hTj ' (1^) 

^± ~ „9ff4-J/.. /.M/9/A ^ N ' \^^) 



e 



2H+J 



:z,/z,y/^{x^-\^) 



in which Zi are defined as the fugacity of species i, that is, Zi = exp(/ij), and A± the 
eigenvalues of the transfer matrix. Note that our choice of (free) boundary conditions is 
expressed by the last term of Eq. [HI through the values of x±. 

Our next step is to determine the size distribution of linear chains of assemblies p{N) 
from Eq. [2] in terms of the concentrations of the two species involved. To do so, we need 
to eliminate the chemical potentials /ii and /i2 from Eq. [2] and Eq. [121 The values of the 
chemical potentials can be established from the conservation of mass for either of the two 
species. 



^ Arp(Ar)z/ = $1 + $2 = $, (15) 

7V=1 

oo 

Y^Np{N)u{S{N)) = <^i-<l>2, (16) 

7V=1 

where $i and $2 sue the molar fractions of species 1 and 2, respectively, and (S) is the 
average spin value, which can be calculated from the partition function as {S{N)) = 
{l/N)dlnZN/dH. Here, we have for convenience introduced the overall molar fraction 
of both species $. Note that the low density approximation used in Eq. [Umplies that our 
results are only valid as long as $i ^ 1. 

To calculate the above sums, we rephrase the partition function Zjv(/ij, &jj, Oj) in terms 
of intensive and extensive parts for the case A^ > 1. Eq. [5] can be rewritten as, 

p(l)z/ = 2;iexp(ai) + Z2exp(a2), (17) 

p{N > l)v = ZN{zi,bi„ai) = J2 Xie-'A^/Xi (18) 

where Aj = Aj exp{b)y/ziZ2 with b = (6ii + 622+26i2)/4. Obtaining a convergent sum requires 
that Aj < 1, which restricts the possible values of the chemical potentials. The two sums in 
Eq. [15] can now be calculated as a geometrical series and be simplified to 

- 2 — A- 

^ = zi exp(ai) + Z2 exp(a2) + ^ Xie'^A] _ ' (19) 



<l>i-<l>2 = 2;iexp(ai)-2;2exp(a2)+ J^ '''''~ \.(i _' A.)2 [~^'(^~^^) dH ^^'^'^'^dW^' 

1=+,— * * 

(20) 
Equations [19] and [20] should be solved simultaneously to determine Zi and the relevant quan- 
tities of interest such as the fraction of self-assemblies, the average degree of polymerization 
and so on. Therefore, these equations are the central equations that are used throughout 
the rest of this paper, when analyzing the model for different free energy parameters and 
concentrations of species. 

To find the simultaneous solution of these two equations for known concentrations of 
both species, we need to look for solutions that satisfy the conditions < Zi < z* < 1, 
where the z* are the maximum value for the fugacity of species i to be discussed in the 
next section. However, to explore the dependence of phase behavior on the concentration 
of the two species, our strategy is to vary < Zi < z* and calculate the corresponding 
molar fractions. Once the chemical potentials from Eqs. [TH] and [2D] are determined, it is 
straightforward to calculate the relevant quantities of interest, such as the number-averaged 
size of assemblies defined as, 

where the sum in the denominator can be easily calculated as a geometrical series. 



oo 



A 
Y^ p(N)u = zi exp(ai) + Z2 exp(a2) + ^ ^'^'''yJiZIT) ^^^^ 

Likewise, one can obtain the fraction of self-aggregates /, defined as 

f = 1 _ ^(^)'^ = 1 _ ^iexp(ai) + 2:2exp(a2) ,^^. 

Having Set up our model for the linear self-assembly of bidisperse monomers, we have 
all the tools in hand for a calculation of the quantities of interest. To gain insight into the 
phase behavior of the system, we first look at a few limiting cases, where analytical progress 
is possible. Therefore, the next section is devoted to the investigation of the self-assembly 
behavior in the low and high concentration limits, and a thorough numerical investigation 
of our model is postponed to section after that, i.e., in section V. 



IV. LIMITING CASES 

A. The monodisperse case 

It is instructive to explore the behavior of system in the limiting case where one can 
simplify the equations and obtain analytical results. The results obtained in these cases 
shed some light on the general "phase" behavior of the system at hand. Before we go into 
details it is important to notice that we recover the results for the monodisperse case (32| 
from the present theory in the limit Z2 — )■ 0, i.e., in the limit of zero concentration of second 
species. In this limit, iJ — > 00, A+ — )■ e"^+^ and A_ — t- yielding A+ — )■ Ziexp^bn). It is 

9 



important to point out that the second eigenvalue of the transfer matrix becomes identically 
zero due to the vanishing concentration of the second species. This is different from the 
usual ground-state approximation valid in the limit of large N, where the contribution of 
the second eigenvalue to the partition function becomes negligible relative to the first. 
We find that the number density of monomers and assemblies are 

p(l)z/ = 2;iexp(ai), (24) 

p(iV >l)z/ = exp(-6n)AA 

where A+ is the solution of the cubic equation of the form Kj^Kai + exp(6ii)A^(2 — A+)/(l — 
A+)^ = $ with Kai = exp(— oi) the activation constant. This is identical to what was 



already known for monodisperse systems, confirming the consistency of our model [32, 134 
The number- averaged length of assemblies for sufficiently high concentrations $i ^ $^ can 
be written as 

W - ^. (25) 

which can be expressed in terms of the molar fraction of monomers in the simple form of 

(iV) ^ V($i/$t-l)/ir„i. (26) 

Here, ^^ = exp(— fon + ai) is a critical concentration beyond which one can consider the 



system in the assembly-dominated regime [33, |3J|. Provided that ai ^ 1, this leads to a 
very sharp cooperative polymerization transition demarcated by ^t. Eg. ( [26]) . and produces 
the well-known square-root growth law of self-assembled polymers |33| . 

The fraction of materials in assemblies in the low and high concentration regions in that 
case obey 

/ ~ 2i^,i$i/$* if $1 < $t (27) 

/ ~ 1 - $*/$! if $1 > $*. (28) 

This shows that for low concentrations the fraction of assemblies grows linearly with concen- 
tration, with a slope proportional to the activation constant. If oi 3> 1 then Kai <^ 1 and 
/ — !■ for $ < $*. On the other hand, for high concentrations the fraction of self-assemblies 
differs from unity by the amount equal to the ratio of the critical concentration to the total 
concentration of monomers. We refer to a recent review by one of us for a more detailed dis- 



cussion of activated equilibrium polymerization |3J| . Next we consider the limiting behavior 
of bidisperse systems. 



B. The bidisperse case for <I> > <!)* 

If the total concentration of monomers is large enough so that (A^) ^ 1, which is the case 
for $ > $* when the activation constants Oj are sufficiently large, we can apply the ground- 
state approximation for the partition function of the Ising model and, therefore, ignore the 
smaller eigenvalue. As a result, the number average length of tapes takes the simple form of 

{N)n ^ \^. (29) 



10 



Furthermore, this leads to a number of simphfications and allows us to get an insight into 
the high-concentration behavior of self-assembly. Particularly, we can find an analytical 
expression for the fraction of self-assemblies. In this limit, the concentration of monomers 
contributing to the assemblies becomes dominant relative to that of free monomers and 
A+ = A+ exp(b)y/ziZ2 approaches its limiting value, i.e., A_|_ — ?> 1. As before, b denotes the 
average of the binding free energies. Therefore, the concentration of free monomers reaches 
its critical value zjj'exp(ai) -|- Z2 6xp(a2) = $*, which corresponds to the maximum molar 
fraction of free monomers. In some specific cases, this maximum molar fraction of free 
monomers characterizes the transition from minimal assembly to an assembly-dominated 
regime. Again, provided that Oj ^ 1, this leads to a sharp cooperative polymerization 
transition. 

The critical concentration depends, in principle, on the stoichiometric ratio of the two 
components, a = <l>i/<l>25 in addition to the various free energy parameters that describe the 
model. Relatively straightforward algebra gives for the fraction of self-assemblies a universal 
curve for high concentrations and of the simple form 

where $*(q;) denotes the critical concentration that now depends explicitly on the stoichio- 
metric ratio a. We can obtain the molar fraction of free monomers for each of the species, 
hence, the critical concentration as a function of a, by solving the two following equations 
simultaneously, 

A4zl4) = 1, (31) 

(9A , , a -I 



dH'''^-' a + V 

where the second equation results from simplifying $i(-2i, Z2)/^2{zl, ^2) = a in the ground- 
state approximation. 

We were not able to find an exact analytical expression for the free monomer molar 
fractions and critical concentration for arbitrary values of a. However, for the special case 
of J S> 1, we did obtain asymptotic expressions for z^ and Zg- The case J ^ 1 occurs either 
if there is a large asymmetry between the two species, or if the two species do not have a 
large affinity to bond to each other. In this limit, one can simplify the above equations to 
obtain an asymptotic solution for the fugacities and therefore, the density of free monomers. 



-J _ ^, 2e^-^(l + a)-a-^a^ + 4(1 + a)e^J 



n ^ ^i — ^ — ' on , :^.4.7 — — (32) 



and 



(2 + 2e4^(l + a) - 2^a^ + 4e^J{l + a) + a (2 + a - ^/ a^ + Ae^-^ {1 + a) 

H - ^- JT7 T^ — (33) 

^ ^ 2e4^(l + a) ^ ' 

Here, $* = exp(— 6jj + Cj) are critical concentrations already defined for the individual 
species. 



11 



For an arbitrary value of J > 0, we can find expressions only in the limits of small and 
large a, corresponding to the respective limits e~^ <^ 1 and e^ ^ 1 in the Ising model. The 
molar fraction of free monomers of either of species in the high concentration regime are, 

$( ~ a^le^-^ <l>^ ~ (1 - a)^l a < e"^"^, (34) 

${~$t(l--) H^-^le^^ a>e^^, (35) 

which agrees with the small and large a limits of Eqs. I32ll33| as they should. 

An interesting conclusion can be drawn if we compare the ratio of molar fraction of free 
monomers of the two species, i.e., a^ = ${/$2 "with a. 

a^/a ~ e^^$t/^2 « < e"^'^ (36) 

a^/a ~ e-^-'^t/*; a > e^-^ (37) 

These results show that, for both small and large a values, the ratio of molar fraction 
of free monomers is different from a, pointing to fractionation effects even in the "high- 
concentration" regime, that is, concentration very much higher than the critical one. We 
will discuss this issue further in the next section where we investigate numerically the full 
concentration behavior of self-assembly. 

As a result, we obtain the following expressions for the critical concentration of the 
bidisperse system. 

$*-$; + a($te'^ -^) tt < e-^^ (38) 

$*-$: + -me""'^ - $:) a > e^^ 
a 

These functional forms clearly demonstrate that the critical concentration not only depends 
on the critical concentration values of the two species, and their ratio a, but is also strongly 
influenced by the value of coupling constant exp(4J). Notice that in the limits a — )■ and 
a — )■ oo, we recover the critical concentrations of species 2 and 1, respectively, verifying the 
self-consistency of our calculation. 

From the above discussion, we recognize four different regimes: a <^ e~^'^ and a ^ e^*^, 
and for both of these, the cases e^-^ < 1 and e^-^ > 1. If a < e~^-^ and e^'^ < '^2/'^i5 the 
critical concentration decreases linearly upon increase of a with a slope proportional to $35 
i.e., $* ^ $2 — a;$2- However, for a <^ e^"^ and e"^-^ 2> $2/*^!' ^^^ slope is determined by 
e^'^, i.e., <l>* ^ $2 + ct^iG^'^ ■ This leads to an initial increase of the critical concentration for 
sufficiently small values of a. 

In the other extreme of a 3> e^'^, we can also make a distinction between the cases 
e^'^ < <l>^/$2 and e^-^ > $*/<l>2- In the case e'^-^ < $1/^25 we expect a decrease of the 
critical concentration for small values of 1/a. On the other hand, if e^"' ^ $]|/$2) the 
critical concentration increases upon a relative increase of the species 2 population. 



C. Bidisperse case for $ < <I>* 

If the total number density of monomers is sufficiently small, only monomers and dimers 
are present in the solution, and this corresponds to investigating the limit Zi = exp(/ij) ^ 1. 
In this case, we can make a Taylor expansion of Eq. [19] and Eq. [20] in terms of fugacities Zi 
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up to the second order, i.e., considering only the contribution of free monomers and dimers. 
This gives 

$1 + $2 ~ Ziexp{ai) + Z2exp(a2) + Aexp{bi2)ziZ2 + 2exp{bu)zf + 2exp(622)-22 5 (39) 

in which the first two terms are the contribution of monomers, while the 3 other terms 
present the contribution of dimers of type 12, 11 and 22. Furthermore we have 

$1 - $2 — zi exp(ai) - Z2 exp(a2) + exp{bii)zl - exp{b22)zl. (40) 

Combining equations [39] and HOI we find that a ~ a-^ . The reason, of course, is that in the 
low concentration regime, the assemblies are mainly in the monomeric regime. 

With these approximations, we can solve for zi and Z2 and calculate the fraction of self- 
assemblies / up to the lowest order in terms of total concentration as well as the ratio of 
the two components a = <l>i/$2, 

/ ^ -— — r^[2exp(6i2 -ai -02) + aexp{bu - 2ai) + a^^ exp(622 - 202)]- (41) 

(1 + a)^ 

Again, we find that the first dominant term is linear in the total concentration. Its slope, 
however, depends on the stoichiometric ratio of the two components as well as on the free- 
energy parameters involved in the system, where / — )■ if the activation energies Oj are large. 
Equation UT] demonstrates that even at the level of dimer formation, the coupling between 
the two components cannot be ignored. As before, we obtain similar results as those for the 
monodisperse case in the limits a — )■ 00 and a — )• 0, as one would expect. If a is large then 
/ should be proportional to $1 and if it is small then it should be proportional to <I>2- 

Now that we understand the behavior of our system in the limits of low and high con- 
centrations, in the next section we examine numerically the full concentration behavior of 
self-assembly in the case of "ferromagnetic" or blocky copolymeric ordering, J > 0. 



V. FRACTION OF SELF-ASSEMBLIES AND CRITICAL CONCENTRATION 

We first focus on the dependence on the total concentration of the fraction of dissolved 
material present in self-assemblies, as we vary the ratio of the concentrations of the two 
components. In Fig. [21 we have plotted the fraction of self-assemblies / as a function of 
the total molar fraction $ for the particular case where we have assumed that 612 = &22- 
Results for two cases are shown. In the first case, the binding and activation free energies 
for the two species are slightly different (corresponding to weak bidispersity) and in the 
second the binding and activation free energies for the two species are considerably different 
(corresponding to strong bidispersity). The chosen values of the activation and binding free 
energies in the first case correspond to oligopeptides consisting of 4 and 5 monomers and 
are sufficiently small to result in a relatively gradual self-assembly transition, while in the 
second case the free energy parameters of the first species chosen are typical values for an 
11-mer [T^. 

First, we discuss the case of weak bidispersity depicted in figures [2^ and [2]d. In this 
case, the fraction of self-assemblies increases gradually from zero to 1 for all ratios of the 
two components, a, as one increases the total concentration of monomers. The bidisperse 
assembly curves are in between those corresponding to the two monodisperse cases. We find 
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FIG. 2: The fraction of self-assemblies / as a function of the overall molar fraction of dissolved 
material $ = $i + $2 at different stoichiometric ratios a = $1/^2 of the two components 1 and 
2. Figures Ela andO* show the results for weak bidispersity (611 = 10, 612 = ^22 = 8, ai = 2.3 and 
a2 = 1.6). The squares on the curves in Fig. [2]a signify the fraction of self-assemblies at the critical 
concentration, associated with each a, as discussed in the main text. Figures [2]c and [2]d graphs 
show the case of a strong degree of bidispersity (611 = 20, 612 = 622 = 8 , oi = 6 and 02 = 1.6). 
The squares on each a curve in [2]c present the fraction of self-assemblies at their respective critical 
concentrations. The solid lines show the slopes for each a calculated based on the Eq. Hi] and show 
a good agreement with numerical results at sufficiently low values of <1>. In the insets, we depict 
the slope of / at low concentrations as a function of a according to Eq. HTJ 



that for all the curves a crossover from an increasing slope for relatively low concentrations to 
a decreasing slope for very high concentrations can be observed, demonstrating a transition 
from minimal assembly to self-assembly dominated regime. The activation free energies 
were chosen not to be very large, so we do not observe a sharp polymerization transition 
from to a non-zero value. In Figj2]D, we have shown the fraction of self-assemblies at 
low concentrations. To verify the validity of the low-concentration expansion, for each 
stoichiometry a, we have also plotted the lines corresponding to Eq. Hi] As can be seen, 
in this regime the fraction of assemblies / grow linearly as a function of concentration with 
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the same slope predicted by Eq. HH A careful look at this figure shows that the slope of / 
versus a is non-monotonic. To highlight this, we have plotted the slope of / versus a in the 
inset of Fig. [2]d. The curves demonstrate the highly non-linear effects that are the result of 
mixing different kinds of assembler units. 

Now, we turn to the case of strong bidispersity (large J) as depicted in Figs. ^ and [21i. 
Here, we also observe that the self-assembly curves for / are between the curves of the pure 
species as we vary the ratio of the two components. However, the behavior of the mixed 
material as shown in the curves for the intermediate concentrations is dissimilar to that 
of the monodisperse solutions. We notice that the fraction of self-assemblies grows at low 
concentrations with a slope that depends on a. Importantly, it grows linearly at very low $ 
and in accord with our estimates for $ <^ $* in Eq. UH as shown in Fig. [2ji. 

However, after an initial relatively steep growth stage, an intermediate stage emerges 
for which the fraction of self-assemblies shows little variation with $. Finally, for large 
enough concentrations the fraction of self-assemblies enters into a third regime of growth 
that follows the behavior predicted for high concentration region $ ^ $*, i.e., deviating 
from unity inversely proportional to $. 

As discussed in the previous section, in an assembly consisting of a single species of type 
i, the transition from mainly monomeric regime to an assembly-predominated regime can be 
demarcated by a critical concentration $* = exp(— 6jj + Oj), which is equal to the maximum 
molar fraction of free monomers reached in the high-concentration regime. Indeed, for large 
enough activation free energies Oj, this corresponds to a sharp transition where the critical 
concentration identifies a transition from no assembly to an assembly-dominated state and 
/($*) ^ 0. For small values of the activation free energy, where the transition is not sharp, 
the critical concentration signifies the crossover from the low concentration regime to the 
high concentration regime, and /($*) has a non-zero value close to 0.5 represented by the 
squares in Fig. |2^, for the / curves of pure assemblies of type 1 and 2. It would be interesting 
to see if this critical concentration can also characterize the transition from minimal assembly 
to self-assembly dominated region in the mixture of two species and, if so, how it depends 
on the stoichiometry a. 

In Fig. |2^ and[2t, we have marked on each curve the points that correspond to the critical 
concentrations calculated according to Eqs. EU We notice that in the weakly bidisperse 
curve, the value of / at the critical concentration is around 0.5, therefore, one can think of 
$* as demarcating the polymerization transition. However, in the strongly bidisperse case, 
the value of / at $* can have any value depending on a. Particularly for large values of 
a, we find / to be close to 1. In these cases, the critical concentration seems to mark the 
onset of the transition from the plateau region in / curves to the high-concentration regime, 
where / ~ 1 - $*/$. 

Therefore, it would be interesting to obtain the full functional dependence of $* on a 
and compare it to the results obtained in the limiting cases of J ::^ 1, and a ^ e~^'^ and 
a ^ e^'^ for any J. Here, we extract the full a-dependence of $* by solving the Eqs. ED 
numerically for two sets of values of free energy parameters corresponding to small and large 
J values as presented in Fig. ^. (for J = 0.5 fcsT) and Fig. 12b (for J = 3 ksT). For the 
case of small J in Fig. |3^, we have shown the results of our estimates of $* for a <^ e"^"^ 
and a 3> e^'^ . We find very good agreement with Eqs. [38] for sufficiently small and large a 
values. 

In Fig. |3]d, on the other hand, we have depicted the full functional dependence of the 
critical concentration on a obtained in the limit of large J based on Eqs. I32ll33l Again we 
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find very good agreement. The general trend that we find is that the critical concentration 
value agrees with that of species 2 at small ratios a = <l>i/$2, and approaches the value of 
the critical concentration of species 1 for large enough a values, as it should be. For the 
small J case, the crossover is not quite monotonic in a, as is clear from Fig. [3^. The reason 
is that if a ^ 1 and e^'^ > $2/*^*) ^^ is the case here, we initially observe a slight increase 
of critical concentration as discussed in the previous section. 

The shown curves demonstrate that the width of the region where the critical concentra- 
tion deviates from that of either of the two species is a strong function of J, and the larger 
the asymmetry between the two species is, the wider is this region. More interestingly, 
for the case of strong asymmetry the critical value is identical to that of species 2 up to 
a ~ 10^. This probably means that in this case the mixing of the two species in assemblies 
is not preferred in some intermediate regions, as we will discuss in a following section. 

We can get an estimate of the width of the critical concentration curve from the limiting 
formulas of Eq. [SSlfor the critical concentration. We can expect the deviation of the critical 
concentration from that of species 1 or 2, when the contributions of the terms proportional to 
$^ and $2 become equal. This gives us two onsets for low and high a values, ai = $2/'^i6~^'^ 
and an = ^^^le'^'^ , leading to a width 2$2/$i sinh(4J). Inserting the values of the free 
energy parameters for the curves shown in Fig. [3l we find good agreement for the estimated 
width (results not shown). 
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FIG. 3: The critical molar fraction ^* as a function of ratio of the two components a = <I>i/<^2 
plotted for a) weak bidispersity: 611 = 10, 612 = 622 = 8 , oi = 2.3 and 02 = 1.6, corresponding 
to J = 0.5 , ^l = 4.528 X 10"'' and <i>2 = 0.00166 . The lines show the approximate functions 
valid for very low and very high a values according to Eq. ^E\ b) strong bidispersity: 611 = 20, 
bi2 = &22 = 8 , oi = 6 and 02 = 1.6 corresponding to a large value of J = 3, <I>^ = 8.315 x 10~^ 
and $2 = 0.00166. The line shows the analytical results obtained for the critical concentration in 
the large J limit based on Eqs. [32]|331 



In conclusion, both the fraction of assemblies and the critical concentration not only 
depend on the free energetic parameters and total concentration, but are also sensitive 
functions of the relative abundance of the two components. This implies that even mild 
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contamination with a chemically distinct species potentially has a large effect on the degree 
of polymerization. In particular, having a critical concentration that depends on a, the 
interesting question that arises is how the fraction of the two components in the assemblies 
differs from a. This is the subject of the following section. 



VI. DISTRIBUTION OF MONOMERS IN THE ASSEMBLIES 

Having a clear picture of the a-dependence of the critical concentration, $*, and the 
fraction of material in assemblies, /, our aim is to get a deeper insight into the composition 
of the self-assemblies that are formed. We would like to know how the two species are 
distributed in the assembhes, and whether the relative population of the two species in each 
assembly is the same as that in bulk solution, i.e., a. Of particular interest is whether or not 
each assembly is formed of one type of species or if both species contribute to the formation 
of each assembly. 

First, let us see if the ratio of the densities of the free monomers is conserved, as we 
increase the total concentration of monomers and assemblies begin to form. In Fig. H] we 
have depicted the relative abundance of free monomers as a function of total concentration 
for different values of a. We find that for very low concentrations the density ratio of free 
monomers is the same as a. However, for larger $ values the density ratio of free monomers 
drops and is considerably lower than a. This implies that a larger fraction of species of type 
1 monomers (the species with greater tendency to self-assemble) contribute to the formation 
of self-assemblies. 

The ratio of free and bound monomers keeps on decreasing with increasing concentration, 
until the total population of free monomers saturates, i.e., when $ 3> $*. Indeed, calculating 
the density of free monomers pf in the large a limit, we find that p{/ P2 = $i/$2 exp(— 4J). 
For large as, we observe an initial increase of the ratio of bound and free monomers in 
the vicinity of $*, followed by a subsequent decrease, i.e., the dependence on a is non- 
monotonic. We note that a-^ = {^{)/{^2) "^ill grow if $2 decreases. For a > 1, especially if 
a is considerably larger than one, there is a lack of monomers of the second species in the 
system and if $ ^ $^, the second species also contributes to the assemblies made mainly of 
the first kind. This effect is stronger if the difference in binding energies, and hence, J, is 
small. Therefore, in such a situation ($2)'^ can become small and lead to an increase of a^ 
around $ fsi $^. 

We can also determine the average fraction of monomers of each species along assemblies 
of arbitrary length A^. The fraction of monomers of species j in a specific assembly of length 
A^ is defined as 9j = (l/N) X]j=i ^i ■ Note that 9j in an assembly of length A^ varies from 
one assembly to another according to Pjq{9j). However, we can calculate its average value 
{Oj)n as a function of A^ from the partition function Zn, i.e., {Oj)n = {l/N)d\n.ZM /d\n.Zj. 
The constraint X]i=i('^i + ^i ) — ^ ^^ each assembly of arbitrary length implies that 
{92)n = 1 — {^i)n, therefore, from now on we focus on {9i)n- 

In Fig. El we have plotted (^i)Ar for the special case of $1 = $2, for both small and 
large J values at several concentrations. This case is particularly illustrative, as deviations 
of the fraction of the two species from 0.5 reflects the deviation of the distribution of the 
species from the original distribution, i.e., $1 = $2- As Fig. [5] indicates, for both small and 
large J cases, the fraction of the species 1, and hence that of species 2, differs from one half. 
This figure clearly demonstrates that short assemblies are made mainly of species 2, i.e., the 
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FIG. 4: a) Ratio of the concentrations of the free monomer of species 1 and 2 a-^ divided by the 
ratio of total density of monomers present in the solution a, as a function of total concentration 
shown for different values of a = $i/$2 • The corresponding a values are depicted in the legends, 
a) Weak bidispersity, with an equivalent J = 0.5. b) strong bidispersity corresponding to J = 3 . 
The dotted lines correspond to concentrations <I> = <^J and <i> = $21 ^s indicated in the figure. The 
free energy parameters used here are the same as those of figures [2] and [3l 



one with less tendency to self-assemble while the dominant population of long assemblies 
are species 1 (the species with a greater tendency with self-assembly). The behavior in the 
case of large J is particularly interesting. We notice that at each concentration, assemblies 
shorter than N^i are made purely of type 2 species and very long assemblies are made 
mostly of species type 1, and that the transition from compositionally pure assemblies of 
type 2 to compositionally pure assemblies of type 1 is sharp. This clearly demonstrates the 
"demixing" effects occurring due to the large asymmetry of the two species. 

The behavior that we find is consistent with that obtained, in a more general sense, 
from the 1-D Ising model. This more general picture arising from the Ising model is that 
for very short chains, that is, compared to the correlation length ^0 = exp(2J)/2, there is 
less combinatorial entropy available, simply because there is not enough room to move the 
domains about. As a result, states of either spin up or spin down (in our case assemblies 
merely consisting of type 1 or 2) are preferred. 

In the other extreme of chain lengths much longer than the correlation length ,^0) the 
combinatorial factor can benefit from a large number of fragmentations, leading to many 
bound domains of moderate lengths of the order of the correlation length. For J = 0.5/cb^, 
^0 — 1-4 while for J = Sfc^T, ^0 — 202. Therefore, for small J, almost for any assembly 
length we are already in the regime that mixed assemblies are favored, while for large J, 
except for very long assemblies and sufficiently high concentrations, we are in the region 
A^ < ^0) therefore, pure assemblies are preferred. 

Having discussed the average fraction of the two species in the assemblies {9i)]y, we 
next consider the number-averaged size of the assemblies, (N), as a function of the total 
concentration of monomers. See Fig. |6l In this figure we also show the concentration 
dependence of (A^) of the monodisperse species for comparison. As expected, the average 
length of assemblies increases as we increase the total concentration of monomers. We find 
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FIG. 5: The average fraction of monomers of type 1 for a composition of a = 1 as a function of 
N shown for different concentrations. The upper graphs shows a weakly bidisperse case, with an 
equivalent J = 0.5 , while the lower graph shows a strongly bidisperse case, corresponding to J = 3 
. The free energy parameters used here are the same as those of figures [2] and [3j 



that the average length of the assembUes in the bidisperse case lies in between that of the 
pure species. This implies that, e.g., for a = 1, the mean degree of polymerization obeys 
an apparent growth law that deviates from the usual square-root law over, say, a decade 
in concentration. While for small J, (N) grows monotonically with concentration, for the 
large- J case (N) remains fairly constant at intermediate concentrations due to the demixing 
effect discussed earlier. This means that over a range of concentrations the mean aggregation 
number does not grow at all. 

In this range of concentrations, the species 1 monomers have already formed relatively 
long assemblies, while the species 2 are still mainly in the form of free monomers. The 
average aggregation number, (A^), reflects the average number of assemblies formed by 
both species. Therefore, in the averaging, the growing number of long assemblies made 
of species 1 are taken with the growing number of very short assemblies of the species 
2, and effectively (A^) does not grow. This behavior is seen even when the relative 
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population of second species is small a = 20, this implies that impurities can strongly 
affect the growth of linear supramolecular assemblies. It shows again the important effect 
that impurities can have on observed quantities such as the average degree of polymerization. 
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FIG. 6: The number-averaged degree of polymerization of linear assemblies as a function of 
the overall concentration for different stoichiometric ratios a. a = and a = oo correspond to 
monodisperse cases of type 1 and 2, respectively, a) "Weakly bidisperse case, with an equivalent 
J = 0.5 . b) Strongly bidisperse case corresponding to J = 3 . The free energy parameters used 
here are the same as those of Fig. [2] and [3l The vertical dotted lines in each figure show the 
concentrations $ = ^^ and <I> = <I'2, respectively. 



VII. CONCLUDING REMARKS AND OUTLOOK 

To summarize, we have set up a model that allows us to investigate theoretically the effects 
of molecular bidispersity in dilute solutions of quasi-linear, self-assembling objects. Mapping 
our model onto the 1-D Ising model, we find that different arrangements of the two species 
in the assemblies take place, depending on the relative values of the binding free energies 
involved in the binding. These morphologies correspond to ferromagnetic ("blocky"), anti- 
ferromagnetic ( "alternating" ) and paramagnetic-like ( "random" ) ordering of the two species 
in the polymeric assemblies. See Fig. 1. 

Analyzing our model for the case of ferromagnetic-like (or blocky-type copolymeric) or- 
dering, we find a range of interesting phenomena. The fraction of self-assemblies and the 
value of the critical concentration, which quantifies the crossover to the high-concentration 
regime, not only depend on the free energy parameters of the two components in the system 
but also on the relative abundance of the two species. The degree of asymmetry of the 
two species, described by a coupling constant J, strongly influences the dependence of the 
critical concentration on the density ratio of the two components. 

The larger the value of the coupling parameter J, the wider the region over which the 
critical concentration is different from either that of the two species. Looking at the distri- 
bution of the monomers of each species, we find that a large asymmetry encoded by a large 
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value of J gives rise to a larger asymmetry in the composition of the assemblies formed. For 
sufficiently large values of J (signifying a large degree of molecular asymmetry) this leads 
to the emergence of a demixing region, where pure assemblies, made up almost entirely of 
either of the two species, coexist. 

It is of interest, we reckon, to discuss the analogy of the effective coupling constant J 
in our bidisperse model, and the Flory-Huggins binary interaction parameter x, appearing 



in polymer solutions and mixtures |39|. In a polymer solution, the interaction parameter x 
determines the miscibility of the solvent and the solute. Apparently, the parameters J and x 
govern ordering processes in mixtures, in the former that in supramolecular polymers and in 
the latter polymer solutions or blends. In our model, increasing J leads to the segregation of 
species into compositionally pure assemblies, which is analogous to having a large x between 



binary polymer mixtures leading to macroscopic phase separation [36 . 

It is sensible comparing the results of our study with that of the existing phenomenology 
of polydispersity effects on the phase behavior of thermodynamic systems. Commonly, the 
introduction of polydispersity causes a range of new features in the phase behavior. Two 
important observed effects of polydispersity are: i) A strong widening of the coexistence 
region, i.e., of the density range within which two or more phases coexist; and ii) frac- 
tionation, meaning that the coexisting phases have different concentrations to the different 
particle species present J40|. 

Here, for the crossover from a mainly monomeric regime to a self-assembly dominated 
regime, characterized by a critical concentration, we find a strong dependence of the critical 
concentration on the ratio of the two concentrations. Furthermore, we observe an analogous 
effect to fractionation. The relative composition of the two species in the assemblies and 
hence also that of the free, unbound monomers in solution is different from the original 
parent ratio of monomers put into the system, especially for concentrations larger than the 
critical concentration. 

Finally, although we have investigated the effect of polydispersity for the simplified case 
of bidisperse solutions, we believe that the insights obtained provide a general insight into 
the general phenomenology and the qualitative features of more general assembly behavior 
in self-assembling systems. In the general polydisperse case, one would expect the value of 
the critical concentration, which demarcates the transition from the monomeric regime to 
that where self-assemblies predominate, to depend on the exact form of the distribution of 
the polydispersity attribute, and to be strongly affected by the width of the distribution. 

One may deduce that for a weakly polydisperse system, i.e., one with a narrow distribution 
of the pertinent polydisperse attribute, a monotonic crossover from the monomeric regime to 
the self-assembly regime can be expected. For a strongly polydisperse system, characterized 
by a wide distribution function, the appearance of multiple regions where pure assemblies 
of single species coexist should be expected. Moreover, from our calculations we conclude 
that independent of the degree of polydispersity, monomers with lower tendency to self- 
assemble are found more abundantly in shorter assemblies, while those with a larger affinity 
to self-assembly have a greater contribution to longer assemblies. 

A model to investigate the generalized case of a continuous distribution of polydisperse 
attributes, similar to that of the polydisperse lattice gas model for the liquid-vapor phase 
equilibria, [41|, and its consequences on self-assembly is, currently under development J42[. 
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List of symbols 

bij-. free energy of the bonded interaction between two monomers of type i = 1,2 and 
j = 1,2, in units of the thermal energy /c^T; 

Oj: activation free energy of species i = 1,2 in units of ksT, with Kai defining the 
activation constant; 

nl the occupation number of species of type j = 1, 2 at site i of a specific assembly; 

Q: grand potential energy in units of ksT; 

Z^: the partition function of an assembly of length iV; 

p{N): the number density of assemblies of degree of polymerization A^; 

u: interaction volume; 

(A^) : the number-averaged degree of polymerization N; 

J = |(6ii + 622 — 2612): The effective coupling constant in the Ising mapping; 

H = |[(&ii — 622) + ifJ'i — fJ'2)]'- the magnetic field in the Ising model; 

b = I (611 + 622 + 2612): average binding free energy; 

/ij: chemical potential of species i = 1, 2 ; 

Zi = exp(/ij): The corresponding fugacities of species i; 

X±: the eigenvalues of the transfer matrix of Ising model; 

A± = X±exp{b)^/zlZ2^. The effective fugacities of the bidisperse system; 

<I>j: molar fraction of molecules of species i = 1,2; 

$ = $1 + $2: the overall molar fraction; 

a = <I>i/<I>2- the ratio of molar fraction of the two species ; 

/: mean fraction of self-assemblies; 

$* = exp(— 6j + tti): critical molar fraction associated with species i, demarcating the 
transition from minimal assembly to assembly-predominated regime; 

^*{a): critical molar fraction of the bidisperse system, whose value depends on the 
ratio of molar fraction of the two components; 

${: the molar fraction of free monomers of species i ; 

p-f = ${/$2: the ratio of the molar fraction of free monomers; 

9j = l/A^X]j=i^i ■ the fraction of species of type j = 1,2 in a specific assembly of 
length N ; 
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{9j)n'- the average fraction of species of type j = 1, 2 in the collection of assemblies of 
length A^; 
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